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Abstract 

A mathematical model is proposed to predict the transfer mechanism of calcium salts through a nanofil- 
tration membrane. The model is a combination of Nemst-Planck and film theory equations and it is character¬ 
ized by three transport parameters: solute permeability P s , reflection coefficient a and film thickness 8. In the 
present work, the influence of including concentration polarization phenomena on salt rejection is studied and 
the thickness of the boundary layer formed at the feed side adjacent to the membrane, 8, is estimated. Model re¬ 
sults were shown to be in good agreement with experimental data at various operating conditions of pressure, 
temperature and initial concentration. The results show that including concentration, polarization has not a great 
effect on salt rejection for the present operating conditions. Estimated thickness of boundary layer formed near 
the membrane surface was too small and ranged between 10” 13 and 1(T 12 m. In later work, we plan to perform 
some experiments at different operating conditions as higher initial concentration of calcium salts or lower tan¬ 
gential flow to bring the system in polarization concentration and then to apply, verify and validate the model 
proposed. 
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1. Introduction 

One of the innovative technologies that are 
gaining broad application in potable water treat- 

corresponding author 


ment is the membrane process. These processes 
have been utilized for desalination of seawater 
and salt water for several decades, and since the 
early 80 ’s for the removal of suspended and 
dissolved constituents in surface and ground- 
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water. Efficiency, reliability and relatively low 
treatment costs resulted in mass application of 
membrane filtration systems. 

Nanofiltration (NF) is the most recently de¬ 
veloped pressure-driven membrane process for 
liquid-phase separation [1], It has found wide 
applications in many industries such as sepa¬ 
ration of no-sucrose compounds in the sugar in¬ 
dustry [2], recovery of caustic cleaning-in-place 
solutions in the dairy industry [3], removal of 
heavy metals in soil operations [4], separation 
of surfactant in the detergent industry [5], re¬ 
covery of dye from textile effluent [6], treat¬ 
ment of paper mill process water [7] and remo¬ 
val of colour from industrial waste water [8]. 

The properties of the NF membranes lie be¬ 
tween those of non-porous reverse osmosis 
(RO) membranes (where transport is governed 
by a solution-diffusion mechanism) and porous 
ultrafiltration (UF) membranes where separa¬ 
tion is usually assumed to be due to size exclu¬ 
sion and, in some cases, charge effects [9]. 

Their separation mechanisms involve both 
steric (sieving) effects and electrical (Donnan) 
effects. This combination allows NF mem¬ 
branes to be effective for a range of separations 
of mixtures of small organic solutes (either neu¬ 
tral or charged) and salts [10], Various models 
have been proposed in modelling the rejection 
mechanisms of charged (ionic) species by NF 
membranes. Most of these models were based 
on the Nemst-Planck equation. 

The aim of the present work is to propose a 
model for predicting the transfer mechanism of 
divalent ions through NF membranes by in¬ 
cluding the concentration polarization pheno¬ 
mena and to study its influence on the salt 
rejection. This model is a combination of 
Nemst-Planck and film theory equations and it 
is characterized by three transfer parameters 
which are solute permeability Ps, reflection co¬ 
efficient a and film thickness 8. 

In this paper, mathematical formulation of 
the model is reported. Calcium salt rejections are 
estimated and compared to experimental results 
for a binary system. A comparison between re¬ 


sults with and without concentration polariza¬ 
tion phenomena is made and the thickness of 
the boundary layer on the feed side adjacent to 
the membrane 8, is also estimated. 


2. Theoretical aspect 

Extended Nemst-Planck equation is the fun¬ 
damental relationship governing the transport of 
ionic species through membrane pores. It is now 
well understood that rejection of ions by mem¬ 
brane pores is a manifestation of the equili¬ 
brium partitioning behaviour of the ions at the 
pore-solution interface and the ionic fluxes 
generated through convective, diffusive and 
electrostatic potential driven migration forces 
[ 11 ]. 

By considering the motion of compounds 
unidirectional through the membrane and assum¬ 
ing dilute solutions, the transport equation for a 
specie i can be written as [12]: 




dC, FZ i C i cW 

—- +——- 

dx 9 \T dx 


+ J V C,(1 -<r ; ) (1) 


For a binary system, the permeate flow of each 
specie is: 



dC x 

dx 


+ 


FZ,C, dV 
m dx 


+ J v C x (l-a x ) 


( 2 ) 



+ 


FZ 2 C 2 (W 
91 T dx 


+ J v C 2 (1-ct 2 ) (3) 


By rearranging Eqs. (2) and (3) together with 
electroneutrality conditions: 

Z,Z'C,=0 (4) 

X,Z/7. = 0 (5) 

Eq. (1) becomes: 

J X =-P^- + J V C X (1-* X ) (6) 
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where: 
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P and a are constant characteristics of each salt. 
They depend on the permeability of each ion, 
and on its charge density (further details are 
given in a previous work [12]). 

This differential equation (Eq. (6)) is the so- 
called Kedem & Spiegler equation, in which the 
terms on the right hand side represent transport 
due to diffusion (corresponding to the migration 
of the ion due to the concentration gradient 
effect), and convection (related to ion transport 
by the solvent) respectively. 

Eq. (6) is integrated over the membrane 
thickness (0<x<Ax) using the following bound¬ 
ary conditions: 


c,=c; 

o 

ll 

H 

(8a) 

c,=c; 

at x-Ax 

(8b) 

and J x = 

C* J 

•'J V 

(9) 


The solution is: 


C m \~ 


1- 


exp[(l - ct)J y !P s ] 
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p\ 


1 - <7 


( 10 ) 


f 





m j 


( 12 ) 


and so, Eq. (10) becomes: 

(13) 

1 — O’ 

2. 1. Neglecting concentration polarization 

As the system was operated in a closed loop, 
the concentration of ion 1 at the interface of the 
membrane can be taken equal to that in the feed 
solution, i.e.; C = C 0i (see Fig. 1.). 




1 - er .exp 


(1-oVv 



Fig. 1. Concentration profile in the membrane-solution 
system. 


Eq. (13) becomes: 


C o, = 


( 

1 - cr.exp 

i - . 

i 

1_ 

p 

V 

L V 
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( 14 ) 


and the retention of ion 1 is given by: 


where: P. = — 

Ax 

According to Donnan equations [13], the 
concentration of ion 1 at the feed and permeate 
interfaces can be written as follows: 


C 



z x FN¥ d 

m 


(ii) 


R = 1- 


1 - cr 


1 - cr. exp 


(1 ~<T)J V 


(15) 


We note that for the binary solutions, 
R sait= R ca,ion = R anion • This expression is similar 
to those found in the literature [12]. 


) 
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2.2. Considering concentration polarization 

In this case, predicting the solute rejection 
requires a coupled solution of two transport 
models: the first one describes the transfer mech¬ 
anism in the concentration polarization bound¬ 
ary layer near the membrane at the feed inter¬ 
face based on the film theory, while the second 
describes the transport in the membrane pores 
using the extended Nemst-Planck equation. The 
transport phenomenon is described in Fig. 2. 


Film 



tem. 


Concentration polarisation close to the mem¬ 
brane surface is assumed to occur within a boun¬ 
dary film of thickness 8, which is dependent on 
the mass transfer characteristics of the system 
[ 10 ]. 

According to the film theory [13], the con¬ 
centration polarization is given by: 

J\ =~Di^j L + C l J v (16) 

ax 

Integration of Eq. (16) over the thickness of 
the boundary layer 8 gives: 


C m = C pr (C pi -C^ e xA—J v 


(17) 


Eq. (17) is then equal to Eq. (13), and so: 


Pi 
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The rejection of ion 1 can then be written as: 

1 


R = 1 — 


1 + 


1 - a 



' 5' 

/ 

exp 

, V 'A, 

-exp 
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( 1-44 


D, 


i ) 


(18) 


(19) 


3. Results and discussion 

3.1. Model parameters 

Eqs. (15) and (19) are used for the prediction 
of salt rejection, then, to estimate model para¬ 
meters without and with concentration polari¬ 
zation phenomena respectively. A computer 
program in Fortran Power Station language is 


performed to find the smallest value of error (e) 
while varying the thickness of the boundary 
layer 8 and the two transport parameters a and 
P s over appropriate intervals. Results are shown 
in Tables 1 and 2. 

We note that s is given by: 

£ ’ = s| R,h — R e xp I 

N 1 
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Table 1 


Model parameters without concentration polarization 


Solution Concentration, mmol. L 1 Transport parameters, with polarization concentration 


c 

P s 10 6 m.s' 1 

£ 

Ca(CH 3 COO) 2 0.5 

0.3634 

3.073 

0.0167 

1.0 

0.4089 

0.686 

0.0193 

1.5 

0.4650 

0.985 

0.0123 

Ca(N0 3 ) 2 0.5 

0.2003 

1.061 

0.0185 

1.0 

0.1839 

2.226 

0.0351 

1.5 

0.1673 

2.197 

0.0467 

CaCl 2 0.5 

0.2186 

1.119 

0.0171 

1.0 

0.1474 

0.889 

0.0090 

1.5 

0.1592 

0.919 

0.0186 

Table 2 




Model parameters with concentration polarization 



Solution 

Concentration mmol. L 1 

Transport parameters, without 
polarization concentration 



a 

P s 10 6 8 




m.s 

Ca(CH 3 COO) 2 

0.5 

0.3630 

3.069 0.0163 


1.0 

0.4090 

0.686 0.0192 


1.5 

0.4580 

0.959 0.0123 

Ca(N0 3 ) 2 

0.5 

0.2000 

1.058 0.0184 


1.0 

0.1840 

2.228 0.0350 


1.5 

0.1600 

2.072 0.0323 

CaCl 2 

0.5 

0.2185 

1.119 0.0168 


1.0 

0.1473 

0.888 0.0087 


1.5 

0.1591 

0.918 0.0183 


3.2. Salt rejection with and without including 
concentration polarization 

Estimated model parameters o and P s are 
used to calculate theoretical salt rejection of 
Ca(CH 3 COO) 2 , Ca(N0 3 )2 and CaCl 2 at different 
concentrations. In Figs. 3, 4 and 5, a compari¬ 
son is shown between model predictions and 
experimental results for calcium salts with and 
without concentration polarization. 

Curves show an agreement between model 
and experimental rejection data for all salts as 
the diagram of parity confirms (Fig. 6). 


They also show that including concentration 
polarization phenomena does not have a great 
effect on salts rejection: theoretical and experi¬ 
mental curves are almost superimposed for all 
salts, Fig. 7. This is confirmed by calculated 
model parameter values, with and without con¬ 
centration polarization, which are very close 
(Tables 1 and 2). 

This effect can be explained by both low feed 
salt concentration, which was in the order of 
magnitude of mmol/L, and the effect of per¬ 
mitting tangential cross flow to destroy the con- 
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J v .ltf |m/s| 


Fig. 3. Comparison between calculated and experimen¬ 
tal rejections with and without concentration polariza¬ 
tion phenomena for Ca(N0 3 ) 2 . 



0 5 10 15 20 25 30 


|m/s] 

Fig. 4. Comparison between calculated and experimen¬ 
tal rejections with and without concentration polariza¬ 
tion phenomena for Ca(CH 3 COO) 2 . 

centration polarization boundary layer formed 
at the feed side adjacent to the membrane. 

The estimated boundary layer thickness, 5, is 
too small and ranges between 10 -13 and 1CT 12 m 
for all calcium salts where some examples are 
given in Table 3. 

Nevertheless, we note a slight amelioration 
when including concentration polarization if 
comparing calculated error values reported in 
Tables 1 and 2 respectively. 



|m/s| 


Fig. 5. Comparison between calculated and experimen¬ 
tal rejections with and without concentration polariza¬ 
tion phenomena for CaCl 2 . 



J v .ltf |m/s) 


Fig. 6. Diagram of parity for calcium salts. 

4. Conclusion 

A coupled model, based on Nemst-Planck 
and film theory equations, was proposed to pre¬ 
dict the transfer mechanism of divalent ions 
through nanofiltration membranes. The model is 
characterized by three transport parameters: so¬ 
lute permeability P s , reflection coefficient a and 
film thickness 8. The influence of concentration 
polarization phenomena on calcium salt rejec¬ 
tion was studied. Model results were shown to 
be in good agreement with experimental data 
and show that including concentration polariza- 
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Fig. 7. Theoretical and experimental curves almost su¬ 
perimposed for all salts. 


Table 3 

Estimated boundary layer thickness, 8, at 1.5 mmol/L 


Salt 

Concentration 
mmol. I/ 1 

8. 10 B m 

Ca(N0 3 ) 2 

1.5 

5.5 

Ca(CH 3 COO) 2 

1.5 

9 

CaCl 2 

1.5 

15 


tion has not a great effect on salt rejection for 
the present operating conditions. Estimated 
thickness of boundary layer formed near the 
membrane surface was too small and ranged 
between 10~ 13 and 10~ n m. 

5. Symbols 

C, — Ion concentration, mol. L -1 
C m — Ion concentration at the membrane 
surface, mol. L" 1 

C m _ — Ion concentration in the membrane 

at the feed interface, mol. L _1 
C 0 — Ion concentration in the feed, mol. 

L7 1 


C p — Ion concentration in the permeate, 
mol. L -i 


c; - 

Ion concentration in the membrane 
at the permeate interface, mol. I/ 1 

A - 

Diffusion coefficient of ion i, rrf.s” 1 

F — - 

Faraday constant, C. mol -1 

J, ~ 

Ion i flux, mol. m" 2 .s _1 

A - 

Solvent flux, m.s _1 

p; - 

Ion permeability, m 2 .s~’ 

P s - 

Local solute permeability, m 2 .s“ ! 

9? — 

Universal gas constant, J, KAmol -1 

A - 

Ion i rejection, % 

T 

Absolute temperature, K 

A - 

Ion charge number, dimensionless 


Ion valence number, dimensionless 

x — 

Distance variable, m 

Greek letters 

S — 

Film thickness, m 

- 

Donnan potential, V 

s — 

Error 

a — 

Salt reflection coefficient, dimen¬ 
sionless 


Ion reflection coefficient, dimen¬ 
sionless 

¥ — 

Electrostatic potential of the system. 


V 
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